Coherent destruction of tunneling, dynamic localization and the Landau-Zener 

formula 
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We clarify the internal relationship between the coherent destruction of tunneling (CDT) for 
a two-state model and the dynamic localization (DL) for a one-dimensional tight-binding model, 
under the periodical driving field. The time-evolution of the tight-binding model is reproduced from 
that of the two-state model by a mapping of equation of motion onto a set of SU(2) operators. 
It is shown that DL is effectively an infinitely large dimensional representation of the CDT in the 
SU(2) operators. We also show that both of the CDT and the DL can be interpreted as a result of 
destructive interference in repeated Landau-Zener level-crossings. 
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The coherent control of quantum dynamics of electrons 
by a periodically oscillating external field has been one 
of the subjects of considerable interest both in nanoscale 
solid state physics and in molecular physics under 
laser fields [2]. The interest is now extended to the 
trapped atoms in Bose-Einstein condensates Q, the lo- 
calized spins in molecular magnetsQl, the Cooper pairs 
in Josephson qubits 5], to name only a few. It should be 
noted that, even when the static properties of a quantum 
system is well known, its response to an explicitly time- 
dependent driving field may be nontrivial and, in some 
cases, poses a quite interesting problem. The phenomena 
known as the coherent destruction of tunneling (CDT) @ 
and dynamic localization (DL) Q are such typical non- 
trivial phenomena. Note that the CDT was originally 
found by Grossmann et al. Q for a model of double- well 
potential, but it has been made clear that the essential 
mechanism of the phenomenon can be well understood 
by a two-state model which represents the quantum dy- 
namics between the lowest two states localized to each 
well Bin 

In both CDT and DL, the initial localized quantum 
state never diffuses under a periodic external field. In 
this aspect, these phenomena are similar However, there 
are also some dissimilarities. The DL is an exact result 
obtained in an infinite driven system and is valid irrespec- 
tive of the magnitudes of the transfer matrix element. On 
the other hand, the CDT is derived approximately in an 
extreme case of a small value of the transfer matrix ele- 
ment. In the CDT, the initial distribution is frozen, but 
in the DL, the distribution oscillates around the initial 
value. Thus the relation between the CDT and the DL 
has been controversial [Tl| . 

In this Rapid Communication, we study the relation- 
ship between these remarkable phenomena in a unified 
way. One may study this problem by assuming a tight- 
binding model with a finite length, and by observing the 
change of the behavior of the electron according to the 
change of the chain length[ll|. However, the equation of 



motion for a finite linear chain model does not allow for 
the analytical solution, so that the analysis inevitably be- 
comes a numerical one. We present here a new approach 
to this problem, which shed light upon the internal rela- 
tionship between the CDT and the DL. It will be shown 
that the DL is an infinitely large dimensional represen- 
tation of a generalized version of the CDT, and in fact 
they are closely related each other. 

Let the two-state system |1) and |2) be under an ex- 
ternal field and driven by the Hamiltonian, 



|2>(2|)+7(|1>(2| + |2)(1|), (1) 



where 7 is a constant tunneling matrix element. The 
Schrodinger equation {h — 1), id/ dt\il){t)) — Hi{t)\ip{t)) , 
is cast into the form 

ij^a2{t) = -^a2(t)+7aiW, (2) 

in the representation \ip{t)) = ai{t)\l) +a2{t)\2). Al- 
though this is the simplest equation of quantum dy- 
namics, it cannot be solved analytically for general 
functional forms of E{t). Grossmann and Hanggi 
8] and Llorente and Plata 9] pointed out that 
for a sinusoidal time-dependence of the driving field 
E{t) = EQCOs{ujt), Eq.(I2) is solved approximately in 
the limit of rapid modulation w 3> 7. By substi- 
tuting ai{t) = exp[— i(£'o/2w) sin(a;t)]ci(t), a2(t) — 
exp[i(i?o/2w) sin((jji)]c2(t), Eq.(l2]) is rewritten as 



«^ci(i) = 7exp[i(£'o/w)sin(wt)]c2(t), 

i-^C2{t) = 7exp[-i(i?o/t^) sin(a;t)]ci(t). (3) 

In the limit a; 3> 7, the above equation is integrated 
approximately for a short period I-k jio by assuming that 



2 



ci{t) and C2{t) are constant, since the rapidly oscillating 
terms are separated out as the phase factors. This is the 
inverse adiabatic approximation. We obtain 



il-Ci(T) = -fJo{Eo/uj)c2iT), 
dr 

*^C2(t) = 'jJo{Eo/uj)ci{t), 



(4) 



where r is a coarse-grained time by the unit of 2t:/uj, and 
Jo{Eq/u!) is the zeroth order Bessel function: 
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exp[i(i?o/i^) sm{Luu)]du 



The above equation tells us that the tunneling parameter 
is reduced effectively by the factor Jq{Eq/uj), and even 
vanishes in the case that Eq/lo coincides with a zero of 
the Bessel function. This is the CDtQ. 

An infinite dimensional analogue of the model (1) is 
given by the Hamiltonian 

oc 

H^it) = E{t) ri\n}{n\ 

n— — OO 
oo 

+ A ^ i\n){n + l\ + \n + l){n\). (5) 



n— — oo 



This is a model Hamiltonian for an electron in an infi- 
nite one-dimensional chain under a time-dependent elec- 
tric field, where \n) represents the Wannier state at site 
n. Paradoxically, the Schrodinger equation id/dt\ip{t)) = 
H2(t)\(f{t)) is solved analytically for arbitrary functional 
forms of E(t). We show explicit time-evolution opera- 
tor with a Lie algebra. We define Tp = J2'^=-oo 
T+ = EZ-oo \n+l){n\, T_ = ^^=-00 \n){n + l\ to get 
H2{t) = E{t)To + A (r+ + T_). These operators satisfies 
the relations: 



[ro,T±] = ±r±, [T+,T_] = o. 



(6) 



The solution of the Schrodinger equation is written as 
\ip{t)) = Uit)\^m, where [/(i) = exp+(-i J*H2is)ds). 
By Feynman's disentangling theorem, U (t) is written in 
the form. 



Uit) 



-iA 



(f+(li) +f_(li)) du 



where A{t) 



-iA(u)To 



J* E{u)du, and T±{u) 



T±. Since and r_ 



are commutable, U{t) is rewritten as 

U{t) = cxp[-iA{t)To] exp[-iB{t)] 



(7) 



in which B{t) = A{R{t)T+ + R{t)*T^} with R{t) = 
exp [«j4(u)] du. Since B{t) has the translational sym- 
metry, its eigenstates are given by the plane waves 
|fc) = X^n^*'^"!"') with the time-dependent eigenvalue 



ek{t) = A{i?(i)e-''= -hi?*(t)e''=}. Then the matrix el- 
ement for the transition \n) — > |m) is calculated by using 
the closure relation as 



{m\U{t)\n) 



exp 

Jva- 



—iA{t)m + 1 \ X' 
(2A|i?WI), 



(m - 



(8) 



where x — a'rgi?(t) and Jn(a;) is the nth order Bessel 
function. For a specific choice E(t) — EQCOs{ujt), and 
at each period of the oscillation r — 2ttI/llj {I — 
0,1,2,---), we find A{t) = and R{t) = rJoiEo/uj), 
and the transition probability is given by 



|(m|C/(r)|n)|2 = j2^_„ (2tA| Jo(i?oMI) ■ 



(9) 



This should be compared with the value J^-n (^'''A) 
which corresponds to the case without external field. 
Eq.® indicates that the oscillating external field gener- 
ally reduces the effective transfer by the factor Jo{Eo/uj). 
Especially, if Eq/lo coincides with a zero of Jo (a;), the 
probability to find the electron at site n) oscillates 
temporally and becomes zero, while that to find it at the 
original site n becomes unity at each period 27t/lu. This 
is the dynamic localization (DL)[3l. 

It is clear that the integrability of the Schrodinger 
equation for ([5]) rests upon the commutativity of 
and r_. On the other hand, for the two-state model 
([1]), we can define the analogous operators, Sq = 
i (|1)(1| |2)(2|), S+ = |1)(2|, and S. = |2)(1|, to get 
Hi{t) — i?(i)S'o+7(5'+ + 5'_). These operators, however, 
satisfy a true SU(2) Lie algebra: 



[S'o, S'±] = ±S'±, [5*+, S"-] 2S'o 



(10) 



These are uncommutable relations and the time- 
evolution operator cannot be decomposed in general. 

Let us discuss the relations between the dynamics of 
CDT and DL. We consider the following bosonic repre- 
sentation for ll]) with Schwingcr bosons: 



a At) 



E{t) 



b\h - bib2) + ^ {b\b2 + bihy (11) 



where bi satisfies the commutation relation of indepen- 
dent bosons, = Si j. The Heisenberg equation for 
b\ and &| is given by 



^J^b\it) = -Mb\it)-^bl{t), 
^jblit) = ^bl{t)~^b\it), 



(12) 



which is equivalent to Eq. ^ by the replacement b\ (t) by 
ai{t). The solution of Eq.((T2]) with the initial conditions, 
b\ (0) = b\ , and b^ (0) = b^ is generally written as 



b\{t) 
blit) 



(13) 



Ub = 



a (3 
-(3* a* 
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where a and 13 are time-dependent complex numbers sat- 
isfying |Q!p-|-|/3p = 1. The point is that, if the two-state 
dynamics described by Eq. (2) is solved somehow, it can 
be mapped onto the solution for Eq. (12) and we obtain 
a class of solutions for state vectors in higher dimensional 
representation spaces of SU(2). Recently, Pokrovsky and 
SinitsynfT5| utilized the same argument to derive a class 
of exact formulas describing the nonadiabatic transitions 
for a model of multiple level crossings. 

Let us define the basis states designated by the boson 
numbers, 



l*> = |P,'?> = 



(14) 



where |vac) is the vacuum state of the bosons. The total 
boson number is a constant of motion. We fix p-\-q = 2N, 
and define the site index n hy n = {p — '?)/2. The ba- 
sis states are classified as |n) = \N + n, N — n), {n = 
-N, -N -I- 1, • • • , 0, 1, • • • ,N). The nonzero off-diagonal 
matrix elements are then given by (n + Iji/aln) = 
-fy^{N + n+l){N -n). If we set 7 = A/iV, we have 
a 2N -t- 1-dimensional linear chain model as a represen- 
tation of the SU(2) Hamiltonian, 



N 

H,{t) = E{t) ^ n\n){n\ 

N-1 

+ A /„(|n+l)(n| + |n)(n+l|), (15) 



-N 



in which /„ 
p + q = 2N 



N ) 



Also in the sector 



1 + ^)(1 
1, an analogous expression is obtained. 
Specifically, iov p + q = 1, the two-state model Hi{t) is 
recovered. An important observation here is that, in the 
limit — > 00 with fixed n, the tight-binding model with 
an infinite chain H2{t) is also recovered since /„ — *■ 1. 
Thus the CDT dynamics in Ub can be connected to the 
DL dynamics in the wave functions for the Hamiltonian 

m- 

We now study the time-evolution operator for the wave 
function, V{t), which satisfies, |^'(/)) = V{t)\^) Once 
explicit matrix elements in Ub are obtained, one obtains 
a class of time-evolutions for the driven system (fT5)) . The 
wave function |5'(t)) is given with Ub as 



\^{t)) = ^ (a%\ - f5h\Y (f3%\ + afoty |vac). 

(16) 

By expanding the right hand side, and rearranging the 
terms proportional to b[ 63 , we find the transi- 



tion amplitude for \n 

{m\V{t)\n) = 



{N + m)l{N - my. ,„+„ 
a p 



m ri^m—n 



y {N + ny.{N~ny. 

X E 

r— r„ 

X |a|2(^-'"-'-)( 



N + n 
r 



N -n 
N — n — r 



(17) 



where the summation over r runs from r„i = max{0, n — 
m} to tm ~ min{A^ + n, N — m}. This is rewritten as, 

{m\V{t)\n) 



- y {N + ny{N-ny " 

(18) 

where x = 2|ap — 1, and P^Z^"^~^"ix) is Jacobi's 
polynomial Ul defined as. 



n + a 
n — r 



(x-iYix+iy^^^ 



This expression of V{t) is valid in the region m — n > 
0, m -\- n > 0. In other regions, {m\V{t)\n) is given 
by the replacement; m — > — n, n —>■ — m, a* — > a for 
TO — n>0, m + n < 0, m ^ n, n ^ m, /?*— > — /3 for 
TO — n<0, m + n > 0, and m — > — m, n — > —n a* —> 
a, (3* ^ (3 for m — n < 0, m + n < 0. 

Now set E{t) to a sinusoidal modulation, E{t) = 
Eocos{ujt) with 7 = A/N. The condition for the rapid 
modulation limit uj ^ j is satisfied for Eq. (|12p in the 
limit A^ ^ 1, and it is solved just the same way as the 
corresponding equation for the c- numbers ([2]). We obtain 



a(T) 



exp 



I exp 







t — smujT 


cos ^ 


. 











sin ^ 







1^ 

'A 



Note that the time t is coarse-grained by the unit 2it/u}. 
The following formula is easily proved by using Stirling's 
formula I14I, 



lim N^^P, 



(a,b) 



N 



1 



2N- 



Jai\z\). (20) 



Then, inserting Eq. p^ into Eq. p^ . and noting that 
X — cos {2A/NjQ{EQ/uj)t), we get, in the limit A^ 00, 



(m|y(T)|n) 



exp 



.Eom . TT ■ 

-z smtjr + I — (to — n) 

uj 2 



xJ„_„(2Ar|Jo(£;oMI) 



(21) 



This is exactly the same as the formula ([S]) including the 
phase factor. Especially, when Eq /lu coincides with a zero 
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of Jo{Eo/uj), the CDT occurs in Ub, while the DL occurs 
in V{t). Thus it is shown that the DL is an infinitely 
large dimensional representation of a generalized version 
of the CDT. 

One of the special cases of a class of the Hamiltonian 
([1]) that allows the exact solution is the Landau-Zener 
modelfl^, [l^ E{t) = vt. The solution is written in 
terms of Weber functions, and the transition probability 
from one branch to another according to the temporal 
evolution from t — —oo to t = oo is given by the cele- 
brated Landau-Zcncr formula[T5l[l6|. One of the present 
authors [l7j pointed out that CDT can be regarded as a 
result of destructive interference between the transition 
paths for repeated Landau-Zener level crossings. The 
above result suggests the possibility to extend this view 
to the DL. In the case E{t) = vt {v > 0), the transition 
matrix elements without adiabatic phases are given by 



a 



(22) 



where P = exp [— 27r(5] is the Landau-Zener nonadiabatic 
transition probability with 5 = A'^ /{N'^v), and (j> is the 
Stokes phase given by = 7r/4-|-argr(l — i(5) -|-(5(ln(5 — 1) 
in which T(z) is the F function. The transfer matrix for 
the two-state Landau-Zener model can be mapped onto 
the 2N + 1-site representation S as before. Noting that, 
in the limit iV > 1, ~ 1 - 27r(5 and /? ~ -V^nSe'"'/'^ , 
we find for the matrix element {m\S\n) at the crossing 



(r7i|S'|n) = exp — i — (m — n) 



(2^277 A/ v) . (23) 



This formula agrees with the exact formula obtained from 
Eq.® , as it should. 

For a repeated crossings of the two-state model driven 
by E{t) = Eocos{ujt), and in the case that Eq is much 
larger than 7 and uj, we can approximately decom- 
pose the whole process into sudden transitions at level- 
crossings and the free propagation between them [13] ■ 
The velocity of energy change v is given by the value es- 
timated at the crossings, v = Eqlu. This is also mapped 



onto the 2A'' -I- 1-dimensional representation. Thus, for 
a double crossing within a period of the oscillation, say 
at <i — 7r/2a; and t2 — 37r/2cj, we have the transition 
amplitude {m\T\n) from \n) to |m) in the limit N 00 
as a sum of all contribution from the intermediate states, 



(toITIti) = 



00 

E 



{m\S\l)e-''''{l\S'\n 



(24) 



00 



where S is the transpose of 5', and Q = 
J^^ dtEo cos {ujt) = 2Eq/uj. The summation is carried 
out exactly by using Graf's formula p^. 



E 

m— — 00 



= J, (V^' + C'-2zCcos0) ^ _ 

valid for real numbers z and ^. The transition probability 
is thus obtained as 



z — (e 



2Ay^2LjjTTE^sm 







-1 






UJ 



(25) 



If one notices the asymptotic formula Jo{x) ~ 
■y/2/7ra; sin {x + (7r/4)) for x ^ 1, it can be seen that 
the formula ([9]) agrees at t = 27r/w with the above one 
in the limit Eq/w ^ 1. The phase factor tt/A is noth- 
ing but the Stokes phase at the level-crossing in the di- 
abatic limit. Thus it is revealed that, in the level of the 
two-state model, the CDT is a result of interference be- 
tween the two intermediate transition paths, while in its 
infinitely large dimensional representation, the DL is a 
result of interference between infinite number of inter- 
mediate transition paths. 
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